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The anticommuting properties of fermionic operators, together with the presence of parity conser- 
vation, affect the concept of entanglement in a composite fermionic system. Hence different points 
of view can give rise to different reasonable definitions of separable and entangled states. Here we 
analyze these possibilities and the relationship between the different classes of separable states. We 
illustrate the differences by providing a complete characterization of all the sets defined for systems 
of two fermionic modes. The results are applied to Gibbs states of infinite chains of fermions whose 
interaction corresponds to a XY-Hamiltonian with transverse magnetic field. 



I. INTRODUCTION 



The definition of entanglement in a composite quan- 
tum system [l| depends on a notion of locality which 
is typically assigned to a tensor product structure or 
to commuting sets of observables Various a priori 
different definitions can then be formulated depending 
on requirements concerning preparation, representation, 
observation and application. Fortunately, most of them 
usually coincide. Q 

In the present article we investigate systems of 
fermions where several of these definitions differ due to 
indistinguishability, anti-commutation relations and the 
parity superselection rule. We will provide a system- 
atic study of the different definitions of entanglement 
and determine the relations between them. To this end, 
we will consider fermionic systems in second quantiza- 
tion. That is, the entanglement will be between sets 
of modes or regions in space rather than between par- 
ticles. The latter case was studied in first quantization 
in 0, H, Q whereas entanglement between distinguish- 
able modes of fermions has been calculated for various 
systems in 0, i i [13, El, El ■ 

That the presence of superselection rules affects the 
concept of entanglement has been pointed out and stud- 
ied in detail in [3, [13, [H, El, [i3| ■ There, the existence 
of states was shown which are convex combinations of 
product states but not locally preparable, thus two rea- 
sonable definitions of entanglement already differ. 

In the following the differences will mainly arise from 
an interplay between the parity superselection rule and 
the anti-commutation relation of fermionic operators. 
The different mathematical definitions will carry physi- 
cally motivated meanings corresponding to different abil- 
ities to prepare, use or observe the entanglement, as well 
as to differences between the single copy case and the 
asymptotic regime. 

In section |lT] we introduce the basic ideas and tools 
used in the rest of the paper. We start by defining the 
different sets of product states in section UIIl From them. 



several sets of separable states are constructed by convex 
combination in section HVl It is shown that they all cor- 
respond to four different classes each of which contains 
the previous ones as proper subsets: 

1. States which are preparable by means of local op- 
erations and classical communication (LOCC). 

2. Convex combinations of product states in Fock 
space. 

3. Convex combinations of states for which products 
of locally measurable observables factorize. 

4. States for which all locally measurable correlations 
can as well arise from a state within class 3 above. 

Section |V] analyzes the asymptotic properties of the 
various sets of separable states. As an illustration of all 
these concepts, section IVTl shows the complete character- 
ization of the different sets in the case of a 1 x 1-modes 
system, and their application to the thermal state of an 
infinite chain of fermions interacting with a particular 
Hamiltonian. In order to improve the readability of the 
paper, we have compiled the detailed proofs of all the 
relations in section IVIll 



II. PRELIMINARIES 



The basic objects for describing a fermionic system 
of m modes are the creation and annihilation opera- 
tors, which satisfy canonical anticommutation relations. 
Alternatively, 2m Majorana operators can be defined, 
C2fe_i := a\ -l-afe, C2k ~ Ofc), foi' k = l,. ..2m, 

which satisfy {ci,Cj} = Sij. Each set generates the al- 
gebra C of all observables. A bipartition of the system 
is defined by two subset of modes, A = 1, . . . uia and 
B = niA + 1, . . . m. We will denote by A {B) the op- 
erator subalgebra spanned by the niA (ruB) modes in A 
{B). 
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If rife is the occupation number of the fc-th mode, i.e. 
the expectation value of the operator a\.ak, the Fock basis 
can be defined by 

\n,,...n,^)^ia\r^ ...{alr-\0). (1) 

The Jordan- Wigncr transformation maps the fermionic 
algebra onto Pauli spin operators so that 

C2fe_i = na«aW, C2fe = na«aW. (2) 

1=1 i=l 

The Hilbert space associated to m fermionic modes (Fock 
space) is isomorphic to the m-qubit space. Due to the an- 
ticommutation relations, however, the action of fermionic 
operators in Fock space is non-local. 

For the fermionic systems under consideration, conser- 
vation of the parity of the fermion number, P = Ylk ^k, 
implies that the accessible state space is the direct sum 
of positive (even) and negative (odd) parity eigenspaces. 
Any physical state or observable commutes with the op- 
erator P, so that we can define the set of physical states 

n {p : [p, P] = 0}. 

Correspondingly, Att and will designate the sets of 
local observables, commuting with the local parity oper- 
ators Pa and Pb, respectively. 

We will call an observable even if it commutes with the 
parity operator, whereas an odd observable will be that 
anticommuting with P. Notice that with this nomen- 
clature, odd observables are not the ones supported on 
the odd parity cigcnspace. On the contrary, an observ- 
able with such support will be even in this notation, as 
it commutes with P. It will be convenient to make use 
of the projectors onto the well-defined parity subspaces, 
Pgj.^^. Any state (or operator) commuting with parity 

has a block diagonal structure p = PeP^e + ^oP^o- 
the local subspaces, a parity conserving operator can be 
written A„ = P^A^P^ + F^A^F^. 

One subset of states of particular physical interest is 
that of Gaussian states. They describe the equilibrium 
and excited states of quadratic Hamiltonians. Moreover, 
important variational states (e.g. the BCS state) belong 
to this category. In various respects Gaussian states ex- 
hibit relevant extremality properties [H, [l^ . Fermionic 
Gaussian states are those whose density matrix can be 
written as an exponential of a quadratic form in the 
fermionic operators [20| . 

p = exp (^-^c^Mcj , 

for some real antisymmetric matrix M. The covariance 
matrix of any fermionic state is a real antisymmetric ma- 
trix defined by 

Tfci = ■^i'ciplck, ci]) , 



which necessarily satisfies iF < 1. According to Wick's 
theorem, the covariance matrix determines completely 
all the correlation functions of a Gaussian state. Pure 
fermionic Gaussian states satisfy F^ = —1, and they can 
be written as a tensor product of pur e states involving at 
most one mode of each partition [21| . 



III. PRODUCT STATES 

We start by defining product states of a bipartite 
fermionic system formed by m = mA + rriB modes, where 
ITT- A (it^b) is the number of modes in partition A (B). 

The entanglement of such system can be studied at 
the level of operator subalgebras or in the Fock space 
representation, thus the possibility to define different sets 
of product states. In Fock space, the isomorphism to 
a system of rriA + ms qubits allows separability to be 
studied with respect to the tensor product C^™-* (g)C^™^ . 
At the level of the operator subalgebras, on the other 
hand, one should study the entanglement between A and 
B subalgebras. However the observables in them do not 
commute, in general, and have non-local action in Fock 
space. On the contrary, Att and Stt, i.e. the subalgebras 
of parity conserving operators, commute with each other, 
so that they can be considered local to both parties. It is 
then natural to study the entanglement between them. 



A. General states 



With these considerations, we may give the following 
definitions of a product state. They are summarized in 
Table H 

• We may call a state product if there exists some 
state acting on the Fock space of the form p = 
Pa'^Pb, and producing the same expectation values 
for all local observables. Formally, [22| 

'PO ■= |p : Pb, [PAiB), Pa{B)] = s.t. 

piA^B^) Pa{A^)pb{B^) (3) 
VA, G X, B^ e B^}. 

• Alternatively, product states may be defined as 
those for which the expectation value of products 
of local observables factorizes, 

PI := {p: p{A^B^)^ p{A^)p{B^) 

VA, e X, P. e S.}. (4) 

• At the level of the Fock representation, a product 
state can be defined as that writable as a tensor 
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TABLE I: Relations among sets of product states. 



product, 

V2:={p:p = p^<E>p''}. (5) 

• From the point of view of the subalgebras of ob- 
servables for both partitions, one may ignore the 
commutation with the parity operator and require 
factorization of any product of observables for a 
product state [l^. This yields another set 

rS {p : p{A B) = p{A)p{B) \IAeA,BeB}. (6) 

The two first definitions are equivalent, VQ = VI. 
They correspond to states with a separable projection 
onto the diagonal blocks that preserve parity in each of 
the subsystems. This means that 

a, p—e, o 

is a product in the sense of V2. 

The three remaining sets are strictly different. In par- 
ticular V2 C PI and V3 C "PI, but V3 ^ V2. The 
inclusion 7^2, Vi C Vl is immediate from the definitions. 
The non equality of the sets can be seen by explicit exam- 
ples as those shown in Table HI The difference between 
Vi and V2, however, is limited to non-physical states, 
i.e. those not commuting with parity [23j . 



• ru ■.^vir\n = VQC]ii 

• V2^ ■~r2r\i{ 

• P3^ -.^vinn 

We notice that p E V2t, is equivalent to p = pA® Pb 
where both factors are also parity conserving. 

With the parity restriction, the three sets are related 

by 

= V2^ C VU. (7) 

The proofs of all the relations above are shown in sec- 
tion 



C. Pure states 



For pure states, all Vit, reduce to the same set. If the 
state vector is written in a basis of well-defined parity in 
each subsystem, it is possible to show that the condition 
of VItt requires that such expansion has a single non- 
vanishing coefficient, and thus the state can be written 
as a tensor product also with the definition of V2. 



B. Physical states 

Being parity a conserved quantity in the systems of 
interest, the only physical states will be those commuting 
with P. It makes then sense to restrict the study of 
entanglement to such states. By applying each of the 
above definitions to the physical states, 11, we obtain the 
following sets of physical product states. 



IV. SEPARABLE STATES 



Generally speaking, separable states are those that can 
be written as convex combination of product states. The 
convex hulls of the different sets of product states intro- 
duced in the previous section define then various separa- 
bility sets. Fig.[T]outlines the procedure to obtain each of 
these sets. Table HIl summarizes the different definitions 
and the relations between them. 
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FIG. 1: Scheme of the construction of the different sets. 



General states 



Taking the convex hull of the general product states, 
we define the sets 



As shown in Fig. [Tl we may take the physical states 
that satisfy the definitions for separability introduced in 
the previous subsection, and hence use iSi n 11 as the 
definition of separable states. This yields the sets 

• 5inn EE Slyr, 

• S2'^ ~ 52 n n, 

• 52nn = 52^. 

Only 52',^ is different from the separable sets defined 
above. Actually, given an SI state that commutes with 
P, it is possible to construct a decomposition according 
to Sljr by taking the parity preserving part of each term 
in the original convex combination. Therefore tSl H 11 C 
iSItt, while the converse inclusion is evident. For 53 fl H, 
on the other hand, it was shown in [23j that any parity 
preserving state in 53 has a decomposition in terms of 
only parity preserving terms, and is thus in S3t^. 

All the considerations above leave us with three strictly 
different sets of separable physical states, 



S2-^ (Z 52 ^ G 51-n- 



(8) 



• 51 := CO (PI), 

• 52 co{T2), 

• 53 := co{TZ). 

These contain both physical states, commuting with 
P, and non-physical ones. It can be shown that 53 C 
52 C 51. 

The non-strict inclusion 52 C 51 is immediate from 
the inclusion between product sets. The strict character 
can be seen with an example, in particular in the subset 
of physical states, psi^- 53 C 52 was proved in (23j . 



From the definitions, it is immediate that 527r C 52'7r. 
The inclusion is strict because not every state p G 52'7r 
has a decomposition in terms of products of even states 
(see example ps2\ in Table HI]) . The condition for 52^ 
is then more restrictive. 

From the relation between product sets, 52 C 51, and 
52'7r Q Sl-n. The strict inclusion can be shown by con- 
structing an explicit example of a Vl-n state without pos- 
itive partial transpose (PPT) [23] in the 2x2-modes sys- 
tem. 

The detailed proofs of the equivalences and inclusions 
above are shown in section IVIIBI 



B. Physical states 



C. Equivalence classes 



From the physical sets of product states we define the 
following sets of separable states, 

• 51„ := CO (Plvr), 

• 52^ :=co(P2^), 

Obviously, the corresponding 53^ = 52^. The inclusion 
relations among product states imply 527r C 5l7r. It is 
easy to see with an example that this inclusion is also 
strict. Table HI] summarizes the definitions and mutual 
relations of the various separability sets. 



If one is only interested in the measurable correlations 
of the state, rather than in its properties after further 
evolution or processing, it makes sense to define an equiv- 
alence relation between states by 

pi P2 if pi{A^B^) = p2{A^B^) VA^gX, B^&B^, 

i.e. two states are equivalent if they produce the same 
expectation values for all physical local operators. There- 
fore, two states that are equivalent cannot be distin- 
guished by means of local measurements. 

With the restriction of parity conservation, the states 
that can be locally prepared are of the form 527r, i.e. 
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TABLE II: The different sets of separable states and their relations. 52, contains all states preparable by means of LOCC. 52', 
represents the original definition of separability in Fock space. 51, gives all convex combinations of states for which expectation 
values of products of locally measurable observables factorize. [51,] contains all states which are locally indistinguishable from 
51,. 



P = T,kPk ® Pk^ where [p^^^^\ Pa(B)] = 0. Since the 
only locally accessible observables are local, parity pre- 
serving operators, i.e. quantities of the form /r?(A^_B^), it 
makes sense to say that a given state is separable if it is 
equivalent to a state that can be prepared locally. With 
this definition, the set of separable states is equal to the 
equivalence class of 52, with respect to the equivalence 
relation above. 

Generalizing this concept, we may construct the equiv- 
alence classes for each of the relevant separability sets, 

[5i^] := {p : 3/5 e p}, « = 1, 2', 2. 

From the inclusion relation among the separability sets, 
[52^] C [52'^] C [Sl-^]. And, obviously, C [Si-^]. 

On the other hand, any state p £ [51^] has also an 
equivalent state in S2t^ (see section IVII B|) . so that 

[52,] = [52',] = [51,]. 

This equivalence class includes then all the separability 
sets described in the previous subsection. However, it is 
strictly larger, as can be seen by the explicit example 
P[su] in Table El 



D. Characterization 



It is possible to give a characterization of the previously 
defined separability sets in terms of the usual mathemat- 
ical concept of separability, i.e. with respect to the ten- 
sor product. This allows us to use standard separability 



criteria (see [251 ] for a recent review) in order to decide 
whether a given state is in each of these sets. 

The definition 52', corresponds to the separability in 
the sense of the tensor product, i.e. the standard no- 
tion [H, applied to parity preserving states. 

As convex hull of 7-'2nn, the set 52, consists of states 
with a decomposition in terms of tensor products, with 
the additional restriction that every factor commutes 
with the local version of the parity operator. Using the 
block diagonal structure Pg/jPg-fP^pPQ of any parity pre- 
serving state, each block must have independent decom- 
positions in the sense of the tensor product. Then a state 
will be in 52, iff both P^pP^ and P^pP^ are in 52',. 

A state p is in T'O, if its diagonal blocks are a tensor 
product, 

^ P^®PfpP^®Pf = P^®P^ e7'2,. (9) 

a, P—e. o 

The set 51, is characterized as the convex hull of T'l, = 
7^0,, i.e. it is formed by convex combinations of states 
that can be written as the sum of a parity preserving 
tensor product plus some off-diagonal terms. 

Finally, the equivalence class [51,] = [52,] is com- 
pletely defined in terms of the expectation values of ob- 
servable products ^,_B,. These have no contribution 
from off-diagonal blocks in p, so the class can be charac- 
terized in terms of the diagonal blocks alone. Therefore 
a state is in [51,] iff 

^ P^ ® PfpP^ ® Pf e 52',. (10) 

Of, /3— e, o 
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Since the condition involves only the block diagonal part 
of the state, it is equivalent to the individual separabil- 
ity (with respect to the tensor product) of each of the 
blocks. 



MULTIPLE COPIES 



The definitions introduced in the previous sections ap- 
ply to a single copy of the fermionic state. It is neverthe- 
less interesting to see the stability of the different criteria 
when several copies arc considered, and, in particular, to 
understand their asymptotic behaviour when iV — > 00. 

The criteria 52'^ and 52^ are stable when several 
copies of the state are considered. 



A generic state of a 1 x l-mode system can be written 
in the Fock representation as 



P e 52',, 
p e 52,. 



Moreover, it was shown in [l5| that the entanglement cost 
of 52, converges to that of 52',, so that asymptotically 
both definitions are equivalent. 

On the other hand, 51, and [51,] do not show the 
same stability, although the corresponding individual 
separability is a necessary condition for the separability 
of the multiple copies state. 

• (^su^ pe 5i„ 
. £ ^ P e [51,]. 

It is also possible to prove (see section IVII C|l that 

Therefore, an NPPT state p is also non separable ac- 
cording to the broadest definition [51,] when one takes 
several co pies . This is true, in particular, for distillable 
states [26l.l27l| . This suggests that the differences between 
the various definitions of separability may vanish in the 
asymptotic regime. The strict equivalence of the classes 
in this limit, however, is proved only for the case of 1 x 1 
modes, as detailed in the following section. 



VI. 1x1 MODES 



In the case of a small system of only two modes, it is 
possible to apply all the definitions above to the most 
general density matrix and find the complete character- 
ization of each of the sets. Table IIIII shows this charac- 
terization. 



/l 



V 



X — y + z p q r \ 

p* X — z s t 

q* s* y — z w 

r* t* w* z 



(11) 



where x, y, z are real parameters, and with the additional 
restrictions that ensure p > 0, which include z < x, y, 
and 1 + z > X + y. 

States in VI must satisfy a single relation between 
expectation values, namely (C1C2C3C4) = {C1C2) (c^c^) , 
which reads, in terms of the given parametrization, z = 
X y. This condition is also necessary for states in P2 or 
VS. 

If a state is in 7-'2, it can be written as the tensor 
product of two l-mode matrices, each of them determined 
by one real and one complex parameter. This imposes a 
number of restrictions on the general parameters above, 
that can be read in Table IIIII Since 52 corresponds to 
separability in the isomorphic qubit system, a state will 
be in 52 iff it has PPT 

According to [2^, a state in V3 has zero expectation 
value for all observable products A^B^, and one of the 
restrictions of p to the subsystems is odd with respect 
to the parity transformation. There are then two generic 
forms of a product state V3 in this system, as shown in 
the table. 

If we restrict the study to physical states, i.e. those 
commuting with P, the density matrix has a block diag- 
onal structure, and the most general even 1x1 state can 
be written 



/ 1 — X — y + z 



V 



















(12) 



Particularizing the conditions for general product states 
to this form of the density matrix, where p = q = t = 
w = 0, gives the explicit characterization of the physical 
product states according to each definition. 



In particular, the state (jl2|) is in T'l, iff z = xy. 
Convex combinations of this kind of states will produce 
density matrices that fulfill \s\'^ < z{\ — x — y + z) and 
< [x — z){y — z), and thus have PPT. This shows 
that, for this small system, 51, = 52',. 

The independent separability of both blocks of p, that 
determines separability according to 52,, requires that 
r = s = 0, i.e. that the density matrix is diagonal in this 
basis. 

Finally, the characterization (|10p of [51,] applied 
to (|lip yields the condition that the diagonal of p is sep- 
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TABLE III: Characterization of the sets for a 1 x 1-modes system. 



arable according to the tensor product, so that all even 
states of 1 X l-modes are in [51^]. 

If we look at several copies of such a 1 x 1-modcs sys- 
tem, it is possible to show that 

e [SU] ^ pe S2\. 

Therefore, in this case all the definitions of entanglement 
converge when we look at a large number of copies. 

A. Thermal states of fermionic chains 



All the concepts above can be applied to a particular 
example. We consider a ID chain of N fcrmions subject 
to the Hamiltonian 

F = i ^ (4a«+i + li-c-) - XI 

n n 

+^Y^ (44+1 +h.c.) . 

n 

This Hamiltonian can be obtained as the Jordan- Wigner 
transformation of an XY spin chain with transverse mag- 
netic field [1^, HOl ■ The Hamiltonian can be exactly di- 
agonalized by means of Fourier and Bogoliubov transfor- 
mations, yielding 

, __ N -1 



with Afe = y^[cos2i^-A]'+472sin2^, = 
cosOkUk + isin^fcfl^^,, cos20fe = and = 

We consider the thermal state p = -[^^^^ttj, with in- 
verse temperature /3, and calculate the reduced density 
matrix for two adjacent modes in the limit of an infinite 
chain, by numerical integration of the relevant expecta- 
tion values as a function of the three parameters of this 
model, A, 7 and /3. 

First we may study which values of the parameters 
make the two modes entangled according to each of the 
definitions. As mentioned above, for a 2-mode system 
there is no distinction between the sets Sl-^ and 52'7r. 
Therefore we look for the limits of the separability regions 
52'^ and 52^ for a fixed value of the parameter A. The 
results are shown in Fig. [2l For every value of A we 
may see that the reduced density matrix is in <S2^ only 
if /3 = 0, i.e. for all finite values of the temperature 
two adjacent fermions will be entangled according to this 
criterion. The region S2'j^, on the contrary, changes with 
the parameters, as shown by the plots. 

From a quantitative point of view, the entanglement 
with respect to tS2'^ can be measured by the entangle- 
ment of formation |3l| , 

Ef{p) ^ min Vp,£;(V'i). 

b,'4'i} 

With respect to S2^. it is natural to define the entan- 
glement of formation conforming to parity conservation 
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as 

= min 

U,'4'i} 

where the minimization is performed over ensembles all 
whose tpi have well-defined parity Both quantities 

can be calculated. The results as a function of the tem- 
perature P, for fixed values of A and 7, arc shown in 
Fig. [31 Consistently with the results in Fig. [21 there is 
always non-zero entanglement with respect to 52^, for 
/3 7^ 0. The entanglement of formation with respect to 
iS2'7r is, for any other value of the temperature, strictly 
smaller, and in fact the reduced density matrix starts to 
be entangled at a finite value of /3. 



VII. DETAILED PROOFS 



A.2. 1^2 c Vl 

Proof. The inclusion V2 C T'l is immediate from the fact 
that the products of even observables in the corre- 
spond, via a Jordan- Wigncr transformation, to products 
of local even operators Ae Be in the Fock representation, 
and thus they factorize for any state in 7^2. The strict 
character of the inclusion is shown with an explicit ex- 
ample as pvi; in Tabic [H □ 

A.3. P3 C VI 

Proof. The inclusion Vi C Vl in immediate from the 
definitions of both sets. The example p-pi ^ Vi (Tablc[T| 
shows it is strict. □ 

A.4. V2 7^ 



This section contains the detailed proofs of all the in- 
clusions and equivalences that appear in the text. Ta- 
ble IIVI summarizes all the definitions and the relations 
among sets. 



Product states 



Proof. The example 
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A.l. VQ = Vl 

Proof. States in VQ satisfy the restriction that 

p{A^ B^) = PA{A^)~PB{Brr) 

for some product state p and all parity conserving op- 
erators ^TT, Bt^. Since the only elements or p contribut- 
ing to such expectation values are in the diagonal blocks 
®P|pP^ ®P|, (a, /3 = e, o), the condition is equiv- 
alent to saying that the sum of these blocks is equal to 
the (parity commuting) product state p = pa® Pb- 

The condition for p ^ VI turns out to be equivalent to 
this. We may decompose the state as a sum 



a. P—e,o 



R, 



where p' is a density matrix commuting with Pa and 
Pb, and R contains only the terms that violate parity 
in some subspace. It is easy to check that R gives no 
contribution to expectation values of the form p{ATr i?7r), 
so that p'{At^B.^) = p' (Att) p' (Btt) . On the other hand, 
an operator that is odd under parity has the form 
Ay = FfAyFf+F^AyFf. Therefore p'(,AyBy) = = 
p' {Ay)p' (By). Since p' commutes with parity, and then 
odd observables give zero expectation value, we have 
checked that all expectation values p'{AB) factorize and 
then p' is a product. □ 



fulfills p-p2 G V2, but p-p2 ^ Vi because it has non van- 
ishing expectation value for products of odd operators, 

f-i- (c2C3)pp2 = « 7^ 0. 

On the other hand, it is also possible to construct a 
state as 
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satisfying G V'i (it is easy to check the explicit char- 
acterization for 1 X 1 modes of Table , but p-pj, ^ V2 
because it is not possible to write it as a tensor prod- 
uct. □ 

A.5. V2^ C Vl-^ 



Proof. The non strict inclusion is immediate from the re- 
sult for general states (jA.2p . Actually, the same example 
p-pi is parity preserving and then shows the non eciuiva- 
Icnce of both sets. □ 

A.6. V2^ = V3^ 

Proof. For any physical state [p, P] = 0, the expectation 
value of any odd operator is null. On the other hand, all 
VS states (in particular those in V^^r) fulfill p{AyBy) = 
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FIG. 2: Regions of parameters that correspond to separable reduced density matrices of two neighboring fermions according to 
the different criteria. The different curves correspond to values of the parameters for which one of the conditions of separabihty 
(see Table lln|l is satisfied with equality. For a fixed value of A (A = 0.5 for the left plot, A = 0.95 for the right one), the area at 
the bottom corresponds to values of /3,7 for which the reduced density matrix is in 82'-^. In both cases, the region of 52^ lies 
on the horizontal axis. Notice that for A = 0.95 there is a small range of values of 7 < 0.2 for which the entanglement shows 
up when increasing the temperature, as illustrated quantitatively by the figure below. 
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FIG. 3: Entanglement of formation of the reduced density matrix of two neighboring fermions with respect to the sets cS2, 
(dots) and 52', (crosses), for fixed values of A and 7, as a function of the inverse temperature. 



[23|. Since a state in V2-^ can be written as a product of 
two factors each of them commuting with the local parity 
operator, then the only non-vanishing expectation values 
in these sets of states correspond to products of parity 
conserving local observables. It is then enough to check 
that 

p{A^B^) ^ p{A-„)p{B^) p = pA®PB- 

Given the state p we can look at the Fock representation 
and write it as an expansion in the Pauli operator ba- 
sis, where coefficients correspond to expectation values 
of products cii' eg) ... ® cri™'' . 

Making use of the Jordan- Wigner transformation 



any product of even observables in the Fock space is 
mapped to a product of even operators in the subalgebras 
A, B. So it is easy to see that the property of factoriza- 
tion is equivalent in both languages and thus 



p g p £ Vi 



This equivalence implies also that of the convex hulls, 
52^ EE 53^. □ 



10 



Product states 


Separable states 


Equivalence classes 


General 


Physical 


co(X) 


Physical 


/ u 




50 


= co(PO) 


O U, — O U 1 111 




po, := -po n n 


50, 


= co(PO,) 


/ ± 




51 


= co(Pl) 




[51,] 


PU :=-pinn 


51, 


= co(Pl,) 


P2 




52 


= co(P2) 


52', — 52nn 


[52',] 


^2, :=P2nn 


32-!^ 


= co(P2,) 


52, 


[52,] 






53 


= co(P3) 


53, - 53 n n 


[53,] 


:= P3 n n 


53, 


= co(P3,) 


Relations between sets 


■P2,V3 c ri 

V2^VS 


P2, C PI, 
P2, — P3, 


50 = 51 

52 c 51 

53 C 52 


50, =51, 

*S2 TT G (S Itt 


[50,] = [51,] 
[51,] = [52',] 
[52',] = [52,] 
[52,] = [53,] 



TABLE IV; Summary of the different definitions of product and separable states. 



1. Pure states 



A. 7. For pure states VljT '^==^ 'P2, 

Proof. A pure state |*)(*| G 11 is such that = ±*. 
We consider the even case (the same reasoning applies for 
the odd one). Since such a state vector is a direct sum of 
two components, one of them even with respect to both 
Pa, Pb and the other one odd with respect to both local 
operations, and applying the Schmidt decomposition to 
each of those components, it is always possible to write 
the state as 

I*) =^a,,|e,)|e,)+^A|o,)|0,), 

i i 

where {le.;)} arc mutually orthogonal states with 

pAh) = +\e,) (Psje,) = +|e,)) andjio,)} 

are mutually orthogonal states with PA\oi) — —\oi) 

iPBm = -m). 

The condition of PI, imposes that i?, j^I') = 

(^'jA, 1^1') for all parity preserving observables. 

In particular, wc may consider those of the form 

A. = ^^fc|efc)(efc| +^fc|ofc)(ofc|, 

k 

i?, - ^i?^|£fc)(efc|+B£|0fc)(0fc|. 

k 

On these observables the restriction reads 



Let us assume that the state has more than one term 
in the even-even sector, i.e. ai ^ and a2 ^ (we 
may reorder the sum, if necessary). Then we apply the 
condition to A = A^|ei)(ei|, B = i?||e2)(e2|, and apply- 
ing the equality we deduce jaij^Af |a2|^i3f = 0, and thus 
ai ||a2 1 = 0, so that there can only be a single term in the 
\ei)\ei) sum. An analogous argument shows that also the 
sum of \oi)\6i) must have at most one single contribution, 
for the state to be in Pl^. 

By applying the equality to operators A ~ A°|oi)(oi| 
and B ~ _B5'|£i)(£i| we also rule out the possibility that 

has a contribution from each sector. Then, if G PI,, 
it has one single term in the Schmidt decomposition, and 
therefore it is a product in the sense of P2t^. □ 



B. Separable states 
B.l. 52 C 51 and 52', C 51, 

Proof. The first (non strict) inclusion is immediate from 
the relation between product states IA.2I To sec that 
both sets are not equal, we use again an explicit example. 
It is possible to construct a state in Pl^ C 51^ which 
has non-positive partial transpose and is thus not in 52. 
However, this has to be found in bigger systems than 
the previous counterexamples, as in a 2-mode system the 
conditions for 51^ and 52^ are identical, as shown in 
Table Ell 

By constructing random matrices pa®Pb in the parity 
preserving sector, and adding off-diagonal terms R which 
are also randomly chosen, we find a counterexample 
in a 2 X 2-system such that psi^ G PItt by construction, 
but its partial transposition with respect to the subsys- 
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tern B, p^f , has a negative eigenvalue. 

When taking intersection with the set of physical 
states, the inclusion still holds, and it is again strict, since 
the counterexample psi, is in particular in VIt^. □ 

B.2. 51^ = 5inn 



Proof. Obviously, Slj^ C 51 fl 11. To see the converse 
direction of the inclusion, we consider a state p G 51 nil. 
Then there is a decomposition p = XiPi with pi G VI, 
but not necessarily in 11. We may split the sum into the 
even and odd terms under the parity operator. 



:= E ^4 ^P'^ + PP^P) + E ^4 ^P^ - PP-P^- 



in [iS27r], i.e. that there exists a state in 82^,^ equivalent 
to p. 

For p G [iSItt], there is a p G ^Ijr, i.e. p = ^ ^kPk with 
each pk G T'Itt, producing identical expectation values for 
products of even operators At^B^. If we define 



Pk 



E 

a, /3— e, 



it is evident that p' := ^kp'k produces the same 
expectation values as p for the relevant operators (see 
proof lA.ip . Therefore, p ~ p' . Moreover, since 
p'i,{AyBy) = for all odd-odd products, every pj. G 7^2^, 
and so p' G 527^. □ 



The second term, py, gives no contribution to operators 
that commute with P. Since p is physical, this term also 
gives zero contribution to odd observables, so that 



It only remains to be shown that each pi^ ^{pi+PpiP) 
is still a product state in VItt- But since for parity 
commuting observables all the contributions come from 
the symmetric part of the density matrix, pi^AT^BT^) = 
Pi^iATrBj^), and the condition for VIti holds for pi^. 
Therefore we have found a convex decomposition of p 
in terms of product states all of them conforming to the 
symmetry. 

The analogous relation for 52^^ was shown in [2^. □ 
B.3. 52^ C S2'^ 

Proof. Since 7^2^ = 7'2nn, taking convex hulls and inter- 
secting again with 11 implies that 52^ C S2't^. However, 
not all separable states can be decomposed as a convex 
sum of product states all of them conforming to the par- 
ity symmetry. In particular, the state 



PS2'^ = 7 



/l o\ 

110 

110 

V 1 / 



which has PPT and is thus in 82'^^, is not in S2^ (recall 
that for the 1 x 1-systcm, only density matrices which 
are diagonal in the number basis are in S2t^). □ 

B.4. [51.] = [52'.] = [52.] 

Proof. From the relations 52. C 52'. C 51. and the 
definition of the equivalence classes it is evident that 
[52.] C [52'.] C [51.]. To show the equivalence of all 
sets it is enough to prove that any state p G [51.] is also 



C. Multiple copies 



C.l. 



G 51. ^ p G 51^ 



Proof. An arbitrary state can be decomposed in two 
terms, p = Pe + Po, where 



PE 



and 



Po 



E I 

a, /3— e, o 



E I 

a, 0, 7, 6—e, o 
(a,/3)/(7,«) 



S ■ 



For any state in 51., there exists a decomposition pe ^ 
'Ei^iPE^ PO = J2i>'^Po^ such that p'e + Po e Vl^. 
Let us consider two copies of a state such that p := 
p»2 g gY^^ Then, using the above decomposition of p, 
and taking the partial trace with respect to the second 
system, we obtain a decomposition of the single copy, 
P = PE + Po ^ Ei-^jtr2(p|;) + J2i Aitr2(po)^ Since p^ 
was a tensor product, p|; = p^® Pb^ with A = A1A2, 
B = B1B2, so is tr2(p|,), and therefore p G 51.. □ 



C.2. 



„«)2 



G [51.] ^ p G [51.] 



Proof. Using the same decomposition as above, p = pe + 
po, a state p G [51.] satisfies ps G 52'.. If we consider 
p := p®2 — 4. the condition [51.] on the state of 
the two copies reads 

PE = PE ® PE + Po ® Po ^ 52'., 

in terms of the components of the single copy state. Tak- 
ing the trace with respect to one of the copies, then, and 
using the fact that po is tracclcss, ps G 52'., so that 
pe[5l.]. □ 



C.3. p NPPT 



„«>2 



[51. 
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Proof. Wc may restrict the proof to states such that p G 
[iSItt]- In other case the miphcation foUows immediately 
from the previous resuh (jC.2p . Written in a basis of weU- 
defined local parities, any density matrix that commutes 
with the parity operator has a block structure (analogous 
to that of for the 1x1 case). 



Z' Pee C \ 

Peo D 

poe 

Vet pooj 



(13) 



The diagonal blocks correspond to the projections onto 
simultaneous eigenspaces of both parity operators, p^p — 
(8) Pf pP^ (8) Pf , whereas C = (g) Pf pP^ (E> Pf and 
D = P^ ®PfpP^ «)Pf . 



From the characterization ((10)) of separability, the state 
is in [iSItt] iff all the diagonal blocks pa/s are in 82'^^- 
It is then enough to prove that the partial transpose of 
p is positive iff Pf ® Pf p (g) pFf Pf has PPT. Non 
positivity of the partial transpose of p implies then the 
non separability (iS2'^) of the one of the diagonal blocks 
of p p. 

The partial transposition of the above matrix yields 



r.TB 



( 

p',, C 

(C't) 

\{Dy p'^J 



(14) 



where X' := X^'^ , and the Tb operation acts on each 
block transposing the last uib — 1 indices. 

If we take two copies of the state, we find for the cor- 
responding uppermost diagonal block pee := P^ ® Pf /3 

pPf 0Pf, 



( Pee® Pee C®C \ 

Peo ® Peo D (g) D 

poe ® Poe 

V Ct ® Ct Poo® Poo 



and for the partial transposition 



{Pee) 



Tb 



J 

(15) 



D'®D' \ 

p'eo®P'eo C'®C' 

C'ts5C't p',e®p',e 

\D'^®D'^ p'oo®p'ooJ 

(16) 



/ p'ee®P'ee 



The matrices (|14p and ((T6|) are the direct sum of two 
blocks. Thus they are positive definite iff each such block 
is positive definite. Let us consider one of the blocks of 
P^ . namely 



D' ®D' 



(17) 



Pee®>t 

Let us first assume that p'^^ is non-singular. Applying a 
standard theorem in matrix analysis and making use of 
the fact that our p e [iSItt], so that each diagonal block 
is PPT, we obtain that (I17p is positive iff 

P'ee ® P'ee > ® D'Kp';^^' ® Po";')(^'^ ® D'^), 



which holds iff 



p'ee>D'{pU-'D'l 



Reasoning in the same way for the second block of (I16p . 
one gets that 



[Peef^ >0 4^p^« >0. 



(18) 



The result holds also if the assumption of non-singularity 
of Poo {poe for the second block) is not valid. In that 
case, we may take poo diagonal and then, by positivity 
of (pee)"^^ (or p^'^ for the reverse implication), find that 
D' must have some null columns. This allows us to re- 
duce both matrices to a similar block structure, where 
the reduced poo (Poe) is non-singular. □ 



C.4. For 1 X 1 systems, p®^ G [SU] 



p e S2'^ 



Proof. One of the directions is immediate, and valid for 
an arbitrarily large system, since p G 52'^ implies p®"^ G 
S2't, C ^Itt C [51^]. On the other hand, if we take 
p := p®"^ G [iSItt], then the diagonal blocks of this state 
are separable, in particular P^ O Pf pPf (g) Pf G 52'^, 
which was calculated in p5|) . For the case of 1 x 1 modes, 
with p given by ([T^ . this block reads 








r'\ 





{x - zf 
{s*f {y-zf 



,2 













This is in 52'^ iff it has PPT, and this happens if and 
only if p has PPT, i.e. p G S2\. □ 
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